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Abstract: Chiral gauge groups acting on a lattice fermion field are constructed such that 
all fermion modes (doublers) have the same charge. Details are given for an abelian axial 
gauge group within a perturbative framework. An action based on this group correctly 
reproduces the continuum gauge-current anomaly, while preserving global chiral symme- 
try, locality, rotational symmetry and hermiticity. A Wess-Zumino-like scalar field enters 
naturally to enforce exact chiral gauge invariance. The degeneracy amongst the doublers 
can be lifted as in a vector model. 



1 



The construction of satisfactory chiral lattice gauge models has long been a thorny 
problem. Consequently, there is as yet no lattice formulation of electro-weak interactions, 
even though we've known how to define vector-like models on lattices for a quarter century. 
It was recognized early on [1] and crystallized into theorems [2] that generic models of free 
lattice fermions with certain desirable properties, including chiral symmetry, must have 
equal numbers of right- and left-handed fermions with the same quantum numbers as long 
as the latter are associated with exactly conserved charges. This is the famous doubling 
problem. We then expect according to the theorem that left-right asymmetry can be 
attained if the fermions couple to a non-conserved current. This accords with continuum 
physics. Indeed in a chiral interaction j ■ A, the divergence of the current j will in general 
be afflicted by anomalies. On a lattice, this means that the chiral gauge symmetry must 
be explicitly broken, or equivalently is exact but requires scalar degrees of freedom in 
the action to cancel the dependence on the gauge parameters. The challenge then is to 
carefully break the symmetry so that the anomalies come out correctly without sacrificing 
sacred cows like locality, rotational symmetry and hermiticity. This is what will be done 
here. 

Most approaches to date start with a vector-like model and attempt, with varying 
degrees of success and clarity, to remove or decouple fermionic degrees of freedom to achieve 
left-right asymmetry [3-7]. References [8,9] give critical reviews of these and other methods. 
Perhaps the most promising and scrutinized approach is the Overlap Method [3], inspired 
by the 5-dimensional domain- wall fermion idea [4]. It is fair to say that this and other 
methods so far are technically difficult and not easy to analyze. The very recent approach 
by Liischer [7] proves the existence of abelian lattice chiral gauge models with exact gauge 
invariance when gauge anomalies are cancelled. If true this is a significant achievement. 
Though it is not clear whether this formalism lends itself to practical calculations and how 
easily it extends to non-abelian models. 

The approach given here contrasts with previous ones in that it is based on a chiral 
gauge group, albeit of an unconventional type, and does not depend on the removal of 
fermionic degrees of freedom for its chirality. It treats all doublers identically as in a 
vector model. The doublers thus resemble the families (e, /j,, r) of the standard model. 
One can then deal with the fermion degeneracy the same way as in a vector model like 
QCD: a spin-diagonalization to staggered fermions can be made, or Wilson-like terms can 
be added. It seems possible that the degeneracy can be lifted in such a way that some of 
the doublers may actually be interpreted as families. Smit constructed models pursuant 
to this philosophy [10], but there was no handle on gauge invariance. 

The formalism will be illustrated on an abelian axial-vector model. The gauge group 
and the action for the fermions are constructed below as expansions in the coupling con- 
stant, though it may be possible to find closed form, non-perturbative expressions. The 
gauge current anomaly emerges easily as does the gauge-symmetry-restoring Wess-Zumino- 
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like scalar field familiar from the analysis of anomalous continuum models [11,12]. In pre- 
vious lattice approaches (excepting the Liischer formalism) the gauge symmetry is broken 
more implicitly: the dependence on the scalar field is not easy to extract. 

The main motivation of this work is to lay groundwork for the eventual formulation 
of a lattice gauge model consistent with the phenomenology of electro-weak interactions. 
For this purpose it may not be a drawback if we could not go beyond a perturbative 
analysis in the fermion sector. Non-perturbative electro-weak physics seems confined to 
the Higgs sector. The methodology presented here might lead to insights into how the 
non-perturbative Higgs sector interacts with the fermions. Also, many people suspect that 
the chirality of the SU (2) x U (1) coupling is fundamentally related to the presence of scalar 
fields, the symmetry breaking, and the massiveness of the vector bosons. Perhaps a lattice 
approach, such as the one here, can vindicate this suspicion. 

To simplify expressions I will use the following notation: 



where ip n is a Dirac fermion at lattice site n. Units are chosen whereby the lattice spacing 
a = 1. For ease of exposition I will take the kinetic part of the action to be that of ordinary 
naive fermions: 



As is well known, this action contains 2 d flavors of fermions corresponding to the poles in 
the propagator at the corners of the Brillouin zone. This degeneracy is a consequence of the 
invariance of So under "doubling" transformations, ifj n — > Ttfj n , where T = (— l) n ' i 7 At 7 5 or 
any product of these [13,1]. Including the identity, this is a set of 2 d elements, denoted by 
T^ n \ If ifjn ^ = u(p) exp(ipn) , with p 0, satisfies the continuum Dirac equation, then so 
do the doublers ^ (n) = T^^°\ 

A naive attempt to construct a chiral model, by gauging the axial transformation, 
ip' n = exp(i(775^ n )'0 n , is well known to fail. This is most easily seen from doubling trans- 
formations of the interaction, Si = ^E^ + '( m i m + l^A^im). For half of the T^ n \ Si 
changes sign. Thus half the fermionic modes have charge +g, while the other half have —g. 
It follows that the model is left-right symmetric. This is not what we want. The problem 
can be traced to fact that the naive axial transformation does not commute with all the 

, but only half of them; the rest anticommute. 

One of the central points of this letter is that a unitary, axial gauge group can be 
constructed which commutes with all the and gives rise to a doubler-symmetric action 
with the expected continuum behavior. The simplest infinitesimal transformation with the 
desired property has 8ip n proportional to terms of the form iO^ipn+a, with cr=(±l, ±1, . . .). 
There are 2 choices for the parameters 9 which respect unitarity, namely 9 n + 6 n+a and 




(la) 
(lb) 




(2) 
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@n+a/2- Note that in the second choice, 9 sits on a site of the dual lattice. I will use 
this choice because it gives simpler expressions besides being more elegant. It also reveals 
a kind of duality between vector and axial-vector gauge transformations. (In fact if we 
omit the 75 from the infinitesimal transformation to make it vector-like then half the 
fermions transforms with positive charge and half with negative, just as in the naive axial 
transformation of the previous paragraph.) So to first order the axial transformation is 

V4 = V>n + i(jl5 °n+a/2^n+a (3) 

where I have defined g = g/2 d for convenience because the sum over a has 2 d terms. 
Although this transformation is a matrix in space-time, it commutes with the generators 
of lattice rotations and translations. 

The higher orders in the axial transformation are determined by unitarity and the 
group composition law. In order to form a group it is necessary that the transformation, 
starting at 0(g 2 ), involve the gauge field A^. (This is reminiscent of, though not directly 
related to, Fujikawa's derivation of the anomaly [14], wherein the fermionic functional mea- 
sure becomes dependent on A^.) The group property can also be ensured by introducing a 
group-valued scalar field transforming as r)' n = r\ n + 6 n . However, in the end, the analysis 
turns out to be simpler using just A M , and besides, we want to see how far we can go before 
we have introduce scalar fields. 

Consider 2 successive gauge transformations 

iP' = G[cf>,A]iP 

V>" = G[6, A'W = G[0, A + A(f>]G[<f>, A]ip. 
In order to have an abelian group we must have 

G[d,A + A<f>]G[(f>,A] = G[e + <f>,A]. (4) 

This condition, along with unitarity, uniquely determines the 0(g 2 ) term from the 0(g) 
term. Rewriting the transformation, Eq. (3), with this term now included: 

V'n = Tpn+igi5 Y °k^n+a (5) 

~ \f E Plff {M* - (0 k + 0i)[4m + 5 (01 " O k )]}ip n+p+a + 0(g 3 ) 

where k = n + a/2, / = n + a + p/2 and Ak,i is a sum over As on dual lattice links 
forming a path between k and /. The path may be taken to minimize the deviation from 
the straight-line path between k and /. If there is more than one such path, an average 
is taken. (In Eq. (5) and below, a and p always denote the vectors (±1, ±1, . . .), while p 
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denotes a unit vector.) Unlike the 0(g 2 ) term, the 0(g 3 ) term, which I've calculated, is 
not uniquely determined by unitarity and Eq. (4). There are 2 free parameters. 

To go beyond the perturbative analysis of this paper, it would be necessary to find a 
unitary, closed form solution of Eq. (4), which reduces to the infinitesimal transformation, 
Eq. (3). I have not yet made a serious attempt to do this. In such a pursuit, it may or 
may not be beneficial to recast Eq. (4) as a system of differential equations obtained by 
taking <fi infinitesimal and Taylor expanding. 

It is easy to see that det G = 1 to all orders in g. This is done using det G = 
exp(Tr InG), expanding In G in powers of g, and noting that all odd powers vanish because 
^75 = (as does the trace of the space-time matrix at odd powers). This leaves only 
even powers which being real must also vanish because det G can only be a unit phase 
factor, since G is unitary. Since det G = 1, the fermion functional measure is invariant 
under the chiral transformation, Eq. (5). The anomaly therefore does not appear through 
the variation of the measure: it appears through the variation of the action. 

The local axial gauge group, denoted Uf , defined by Eqs. (3), (4) and unitarity is not 
strictly local. At order g N , G mn has non-vanishing elements out to \n — m\ = 2N in 4- 
dimensions. The fermion action based on Uf will thus have long range interactions at higher 
orders in g. However the 0(g N ) term in the action, tfj m . . .ip n , with max\n — m\ ~ 2N, 
has dimension 4 + N, and is thus suppressed by a factor of a N . We therefore expect no 
violations of locality in processes with energy far below the cutoff, 1/a. 

Now I will construct the action through 0(g 2 ) by applying the gauge principle using 
Ui . Starting from So for naive fermions, Eq. (2), the action, invariant to 0(g) under the 
transformation, Eq. (5), is Sq + Si with 

Si = -ig K l i+) (n,n + fi + a/2)A^n + a/2) (6) 

where the gauge field sits on links of the dual lattice and has the usual transformation 
law: A'^ = A^ + A^9. In the naive continuum limit, we have the usual axial-vector 
interaction. 

The variation of So + Si under the transformation can be written after some algebraic 
reorganization as 

5S 0+ i = \g 2 J2{Kj-\n, n + fi + p + a) [(6 l+fl - 9 k f - ±(0, - 6 k f - \{B l+li - 8 k+ ^) 2 
+ 29 l+ ^(k) - 26 k A^l) + (6 k+f , + 9i +ll )A k+ ^ l+fl - (6 k + 8i)A kyl ] } 

where as in Eq. (5), k = n + a/2 and / = n + a + p/2. There is no S2 depending only on ip 
and A^ which makes 5S vanish. The reason is the anomaly. In order to cancel the (SO) 2 
terms we need 
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Now with S = So+1+2 we are left with the anomaly in the form 

SS = g 2 {K { -Hn,n + v + p + a)[e k J2A + e l+ ^J2 A }} ( 8 ) 

n,fi,p,a Pi P 2 

where Pi (P 2 ) is an ordered, closed path with vertices k, I + fi, I (k,k + p,,l + p). 

In order to expose the anomaly in conventional form, we define the effective action, 
W[A], by exp(-W[A]) = J Dip Dip exp(— S 1 ). It follows that 8 W[A] = W[A+A9] - W[A] = 
— ln(exp(— 5S)) (SS) for small 6>, the average being over the fermion fields. In both 2 
and 4 dimensions, (SS) vanishes identically at 0(g 2 ). This agrees with the fact that in 
2-dimensions there is no gauge anomaly in the axial- vector model. 

To compute (SS) in 4-dimensions in the continuum limit the following easily derived 
ingredients are needed: 

Pi Pi 

(K { ~\n, n + p + a)) = 2igCe IJUX p v (p + a) a F{i v + ... 



where 

C = J (2^ sin2 Pi cos2 Pi cos 4 / ( ^ sin 2 p A ) 2 = 
Then we find the expected result for 16 flavors, 

6W[A] = J cfxOe^F^F^ (9) 

from which it follows that the current 3° — dW/ dA^ has the correct anomalous divergence. 

Here, is, of course, the current coupling to the gauge field, and not the current 
associated with the usual global chiral symmetry Ua, [Vn = ex P('i r Y5 a )^n]- Without mass 
or Wilson terms, remains exactly conserved, keeping the fermions massless. As is well 
known [1], the conservation of owes to the cancellation of the Ua anomalies between 
the doublers. So if we extend this chiral fermion method to the Standard Model, and 
try to interpret some of the doublers as families, then we could get a different pattern of 
Ua anomaly cancellation than is currently understood in the continuum model, where the 
anomalies don't cancel between families. However at this stage this is not an important 
consideration, since experimental signatures for the cancellation pattern, such as baryon 
number non-conservation (and B-L conservation) have not been seen, and are not likely to 
be seen in the near future. 
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To make this axial-vector model exactly gauge-invariant it is not enough to simply 
add other fermions with appropriate charges to cancel the anomaly. This is apparent from 
Eq. (8): (SS) will vanish, but S itself is still gauge- variant. To have exact chiral gauge- 
invariance, a group- valued scalar field r\ n must be introduced transforming as rj' n = r] n + 9 n . 
Then inspection of Eq. (8) shows that 

+k cos[(/(A^?7 - A M )] + 0(g 3 ). 

is gauge-invariant through 0(g 2 ). I added a gauge-invariant kinetic term for r/, which is 
recognized as a fixed-length Higgs term, with r\ being the phase of the Higgs field. In 
view of the analysis above, the linear term in r\ becomes a Wess-Zumino term [15] in the 
effective action W[rj, A^\. 

All mass and degeneracy-lifting terms in this model must be made gauge-invariant with 
the use of the scalar field r}\ neither the mass term ip n ip n nor the Wilson term rt[) n {ip n+ ^ — 
2ifj n + ifj n -^} is gauge-invariant by itself: gauge-invariance requires Yukawa interactions as 
in the Standard Model. It is straightforward, using the chiral gauge transformation, Eq. 
(5), to construct these interactions. For example, the gauge-invariant mass term to 0(g 2 ) 
is 

ij n i(j n - 2ig Y $nl5lpn+*Vn+*/2 ~ 2g 2 ^ ^n^n+a+pVn+a /2Vn+a+p/2 (H) 

which reduces in the naive continuum limit to the expansion of ifj x exp (— 2ig^r] x )^ x to 
0(g 2 ). It is interesting to note that this lattice mass term is not invariant under the U a 
global chiral symmetry, while its naive continuum limit is ( provided rj — > rj + 9 under Ua)- 

Because the action is doubler-symmetric, we can reduce the number of fermion modes 
by spin-diagonalizing in the usual way [16]: a unitary transformation, tf) n — > 7™ 1 . . .74 4 ^ n , 
breaks the action into 4 identical, independent pieces. Tossing 3 pieces away we have an 
action describing 4 Dirac fermions (reduced from 16), with one fermionic degree of freedom 
per lattice site (reduced from 4). Then terms can be added to lift the mass degeneracy of 
these fermions. Interestingly, The interaction Si, linear in A^, in terms of these staggered 
fermions is identical to the linear term in Smit's staggered fermion approach [10] except 
for the coordinates of the gauge fields. Higher order terms are however not similar: Smit's 
action has no gauge symmetry. 

How should we use the action, Eq. (10) ? Anomalous models, such as in Eq. (10) are 
believed to be nonrenormalizable [see, e.g. ref. 12]. To obtain a renormalizable continuum 
model, the gauge anomaly should be cancelled by other fermions, as for example in a model 
with 9 fermion fields, 8 with axial charge, 9=1, and one with q = — 2. To extract physics 
from such a model, the parameters must be near a critical point or surface. Ignoring 
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the fermions for the moment, the Higgs model with just rj and A p has a Higgs phase for 
k > k c with Ma > and a symmetric phase for k < k c with Ma = 0. With an eye 
towards electro-weak physics, let's say we want to obtain a continuum limit in the Higgs 
phase with a 3 parameter model, k, g and y, the coefficient of the mass- Yukawa term 
in Eq. (11). At a point in the Higgs phase near the critical surface, for small g and y 
the truncation of the action at some low order in g should be valid. To calculate some 
observable of ip or one could either perturbatively integrate out the fermions in the 
path integral and numerically integrate (simulate) the resulting expression in A^ and ry, 
or include the fermions in the simulation. One would have to check that answers were not 
sensitive to the order of truncation of the action. 

As we get closer to the critical surface, the fermion mass m/, in units of inverse lattice 
spacing, should approach zero (otherwise there would be no low mass fermion in physical 
units). To get rrif — > we can either tune y to zero as k — > k c or keep y fixed and rely 
on the higher orders in Eq. (11) to reduce rrif. At larger values of y we will need more 
terms in Eq. (11). Similar reasoning might apply to an attempt to remove doublers using 
a Wilson- Yukawa term: the fluctuations in r\ may conspire to send the effective amplitude 
of such a term to zero, giving the doublers a mass much less than the cut-off. This is what 
happens in the Smit-Swift model [17,9]. Unlike in the Smit-Swift model the chirality of 
the model here does not depend on the Wilson- Yukawa term. If need be, the number of 
doublers can be reduced by using staggered fermions or other methods. 

The extension of the above methodology to a V-A abelian model is straightfor- 
ward. We want a transformation, which in the naive continuum limit reduces to tf)'(x) = 
exp[ig9 x ^(l — ^)]ip(x), and which transforms all doublers the same. This is accomplished 
by the following infinitesimal transformation, 

= (! + 7;96n)lfj n - ^<775 J^(#n + 9 n+p )^ n+p . 

p 

The higher orders are constructed using the group property, Eq. (4), and unitarity as was 
done above for the axial transformation. The extension to non-abelian chiral groups is 
slightly more subtle, and will be presented in another publication. 
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